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CLc ABSTRACT 

J-*, , Context: This paper is part of a multi-species survey of line emission from the molecular gas in the circum-nuclear disk (CND) of 

■J-J . the Seyfert 2 galaxy NGC 1068. Unlike in other active galaxies, the intensely star-forming regions in NGC 1068 and the CND can 

~^ ' be resolved with current instrumentation. This makes this galaxy an optimal test-bed to probe the effects of AGN on the molecular 

' medium at ~ 100 pc scales. 

Aims: Single-dish observations have provided evidence that the abundance of silicon monoxide (SiO) in the CND of NGC 1068 is 

enhanced by 3-4 orders of magnitude with respect to the values typically measured in quiescent molecular gas in the Galaxy. We aim 

^ ' at unveiling the mechanism(s) underlying the SiO enhancement. 

^T) ' Methods: We have imaged with the IRAM Plateau de Bure interferometer (PdBI) the emission of the SiO(2~l) (86.8 GHz) and CN(2- 

\>0 ^ 1) (226.8 GHz) lines in NGC 1068 at ~ 150 pc and 60 pc spatial resolution, respectively. We have also obtained complementary IRAM 

r\l . 30m observations of HNCO and methanol (CH3OH) lines. These species are known as tracers of shocks in the Galaxy. 

Results: SiO is detected in a disk of ~400pc size around the AGN. SiO abundances in the CND of ~(l-5)xlO"' are about 1-2 orders 
of magnitude above those measured in the starburst ring. The overall abundance of CN in the CND is high: ~(0.2-l)xl0"'. The 
abundances of SiO and CN are enhanced at the extreme velocities of gas associated with non-circular motions close to the AGN 
(r <70 pc). On average, HNCO/SiO and CHsOH/SiO line ratios in the CND are similar to those measured in prototypical shocked 
regions in our Galaxy. However, the strength and abundance of CN in NGC 1068 can be explained neither by shocks nor by PDR 
chemistry. Abundances measured for CN and SiO, and the correlation of CN/CO and SiO/CO ratios with hard X-ray irradiation, 
suggest that the CND of NGC 1068 has become a giant X-ray dominated region (XDR). 

Conclusions: The extreme properties of molecular gas in the circum-nuclear molecular disk of NGC 1068 result from the interplay 
ji^ ■ between different processes directly linked to nuclear activity. Results presented in this paper highlight, in particular, the footprint of 

$_i( ' shocks and X-ray irradiation on the properties of molecular gas in this Seyfert. Whereas XDR chemistry offers a simple explanation 

C^ I for CN and SiO in NGC 1068, the relevance of shocks deserves further scrutiny. The inclusion of dust grain chemistry would help 

solve the controversy regarding the abundances of other molecular species, like HCN, which are under-predicted by XDR models. 

Key words. Galaxies:individual:NGC 1068 - Galaxies:ISM - Galaxies:kinematics and dynamics - Galaxies:nuclei - 
Galaxies:Seyfert - Radio lines: galaxies 

1 . Introduction strong radiation fields (UV and X-rays) and the injection of me- 

,^^ , .... . . , , • chanical energy (by gas outflows or jet-ISM interactions) are key 

Nuclear activity and intense star formation can shape the excita- j^ g^ients in the feedback of activity. Molecular line observa- 
tion and chemistry of molecular gas in the circum-nuclear disks ^-^^^ ^^^ ^^ ^^^^ ^^ ^^^^ ^^^ ^^^^^^ ^^ ^^^ dominant source of 
(CND) of galaxies. The effect on molecular gas properties of ^^^^^^ -^ g^l^^^^ ^^^^^^ ^^^^ ^^^^^ ^^^^^^-^ ^^^^^. ^^q^^ ^^^ 

y, „ ' T ~ ■ . . -.u .u Tn Aivj T,i » J Ti intense Star formation . Furthermore, nearby galaxies can serve 

Based on observations carried out with the IRAM Plateau de Bure ,, , ^,. ,-,, , 

Interferometer IRAM is supported by INSU/CNRS (France), MPG ^^ local templates of distant galaxies where these phenomena 

(Germany) and IGN (Spain) can be deeply embedded. 
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Some of the tracers of the dense molecular gas phase, which 
are commonly used in extragalactic research, can be heavily af- 
fected by the feedback of AGN activity. In particular the lines of 
HCN (the most widely observed tracer of the dense molecular 
medium) appear as over-luminous in the CND of many Seyfert 
galaxies with respect to other tracers of the dense molecular gas 
(Tacconi et al. 1994; Kohno et al. lOOT: Usero et al. '2004, here- 
after U04; Krips et al. 2008 ). Doubts have been thus cast on 
the universality of the conversion factor between the luminosity 
of the HCN(l-O) line, Lhcn(I-O), and the mass of dense molec- 
ular gas in active galaxies (Gracia-Carpio et al. 120061 120081 
Garcia-Burillo et al. 2006; Krips et al. 2008). X-rays can effi- 
ciently process large column densities of molecular gas around 
AGNs, producing X-ray dominated regions (XDR). It has been 
proposed that the abundances of certain ions, radicals and molec- 
ular species, like HCN, can be enhanced in XDR (Lepp & 
Dalgai-no [T9961 Malo ney et al. |T9961 Meijerink & Spaans 2005; 
Meijerink et al. l2007l l. However, it is still controversial whether 
pure gas-phase XDR models are able to enhance HCN abun- 
dances in AGN to the level imposed by observations. 

The inclusion of dust grain chemistry, not taken into ac- 
count by gas-phase XDR schemes, could solve the HCN con- 
troversy in AGN by enhancing the abundance of this molecule. 
Besides HCN, other molecular species can also undergo signif- 
icant changes in their abundances due to dust grain processing. 
X-rays can evaporate small (~10 A) silicate grains (Voit |199TT ). 
This can increase the Si fraction in the gas phase and then 
considerably enhance the abundance of SiO in X-ray irradiated 
molecular gas (Martin-Pintado et al. 2000; U04; Garcia-Burillo 
et al. 120081 Amo-Baladron et al. 2009). Furthermore, mechan- 
ical sputtering of dust grains in molecular shocks is an addi- 
tional source of dust grain chemistry. Large SiO abundances 
have been found in the nuclei of a number of non-AGN galaxies 
where interferometer maps have been key to identifying large- 
scale molecular shocks (Garcia-Burillo et al. 120001 1200 li Usero 
etal.EOOll). 

To address the role of X-rays and shocks in shaping the 
chemistry of molecular gas in AGN requires a multi-species ap- 
proach. Unlike SiO, other tracers of the dense molecular gas, 
like CN, are not expected to be highly abundant in shocks 
(Mitchell 1 1 9841 Fuente et al. 120051 Rodriguez-Fernandez et 
al. 120 101 ). However, there is a consensus supported by obser- 
vations and theoretical models that the CN radical is a priv- 
ileged tracer of highly ionized molecular gas, typical of pho- 
ton or X-ray dominated regions (PDR or XDR) (Boger & 
Sternberg ESSl Fuente et al. fT993l l2008l Janssen et al. [19951 
Le pp & D algarno l 19961 Meijerink & Spaans l2005l Meijerink et 
al. 2007). Observations of tracers like SiO and CN can then be 
used to quantify the relevance of shocks and XDR chemistry in 
AGN. Furthermore the high-spatial resolution provided by in- 
terferometers is paramount to spatially discriminate between the 
different chemical environments (star formation vs AGN activ- 
ity) that can coexist in a single galaxy. 

NGC 1068 is the strongest nearby Seyfert 2 galaxy, and as 
such is a prime candidate for studying the feeding and the feed- 
back of activity using molecular line observations. Tacconi et 
al. ( 11994 ) and Schinnerer et al. (2000; hereafter SOO) used the 
Plateau de Bure Interferometer (PdBI) to map at high-resolution 
(1^") the emission of molecular gas in the central r ~ 1.5-2 kpc 
disk using HCN and CO lines. The CO maps spatially resolve 
the distribution of molecular gas in the disk, showing a promi- 
nent starburst ring (hereafter SB ring) of ~l-1.5kpc-radius, 
which contributes significantly to the total CO luminosity of 
NGC 1068 (see also Planesas et al. [19911 and Baker ,2000. here- 



after BOO). Furthermore, CO emission is detected in a central 
r ~200pc CND surrounding the AGN. In contrast to CO, emis- 
sion from the HCN(l-O) line, which is also spatially resolved in 
the Tacconi et al.'s maps, arises mainly from the CND. Tacconi 
et al. ( 119941 ) derived a high HCN/CO intensity ratio (~1) in the 
CND. This is about a factor of 5-10 higher than the average ratio 
measured in the SB ring (Usero et al. in prep). Radiative trans- 
fer calculations based on multi-line observations of HCN and 
CO showed that the abundance of HCN relative to CO was sig- 
nificantly enhanced in the CND: HCN/CO~10"^ (Sternberg et 
al. [T994l ). 

U04 used the IRAM 30m telescope to observe with low to 
moderate spatial resolutions (~ 10-30") the emission of eight 
molecular species in the CND of NGC 1068. The global analy- 
sis of the survey, which includes several lines of SiO, CN, HCO, 
H'^CO+, H'2C0+, HOC+, HCN, CS, and CO, indicated that the 
bulk of the molecular gas emission in the CND of NGC 1068 
could be interpreted as coming from a giant XDR. More recently, 
Perez-Beaupuits et al. ( 120071 120091 ) analyzed new single-dish 
data obtained for several rotational lines of HCO^, HNC, CN 
and HCN, and came to similar conclusions as U04. However, 
the low spatial resolution of these single-dish observations are 
not optimal to precisely disentangle the contributions of the star- 
burst and the CND in NGC 1068. 

In this paper, we use the high-spatial resolution (~l-3") af- 
forded by PdBI to map the emission of the (ii=0, 7=2-1) line of 
SiO and the N-2-\ transition of CN in the central r ~ 1.5-2 kpc 
disk of NGC 1068. The spatial resolution of the new observa- 
tions, an order of magnitude higher than that of the 30m sur- 
vey of U04, allow us to neatly separate the emission of the SB 
ring from that of the CND. Furthermore, the SiO and CN PdBI 
maps, together with complementary single-dish data obtained in 
CH3OH and HNCO lines, are used to study the chemical differ- 
entiation inside the CND through an analysis of line ratio maps. 
Line ratios are interpreted with the help of one-phase large ve- 
locity gradient (LVG) models. While the emission of the differ- 
ent lines and species analyzed in this work likely comes from re- 
gions characterized by different physical conditions, the adopted 
one-phase approach allows us to explore the existence of overall 
trends in the chemical abundances of these species, using a mini- 
mum set of free parameters in the fit. We explore the dependence 
of line ratios with the illumination of molecular gas by the X-ray 
AGN source. We also explore the link between shock chemistry 
and gas kinematics in the CND. 

We describe in Sect. |2] the observations, including high- 
resolution SiO, CN and CO maps obtained with the PdBI, single- 
dish data obtained with the 30m telescope as well as X-ray im- 
ages taken by Chandra of NGC 1068. Sect.|3]presents the contin- 
uum maps derived at 86.8 GHz and 226.8 GHz. The distribution 
and kinematics of molecular gas derived from SiO and CN are 
described in Sects. |4] and |5| SiO and CN abundances are dis- 
cussed in Sect. |6l In Sect. |7] we interpret the line ratio maps in 
terms of two different chemical scenarios (shocks and XDR). 
The main conclusions of this work are summarized in Sect. [8] 
We assume a distance to NGC 1068 of D -14 Mpc (Bland- 
Hawthorn et al |1997l ); this implies a spatial scale of ~70 pc/". 



2. Observations 

2.1. Interferometric maps 

Observations of NGC 1068 were carried out with the PdBI 
array (Guilloteau et al. 119921 ) between December 2004 and 
February 2005. We used the BC configurations and six anten- 
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nas in dual frequency mode. We simultaneously observed the 
(v^O, 7=2-1) line of SiO (at 86.847 GHz) and the N^2-l tian- 
sition of CN. The CN(2-1) transition is split up into 18 hy- 
perfine lines blended around three groups (7=5/2-3/2, J-3/2- 
1/2 and 7=3/2-3/2). We have observed the two strongest groups 
of lines at 226.9 GHz (7=5/2-3/2) and 226.7 GHz (7=3/2-1/2), 
hereafter referred to as high frequency (HF) and low frequency 
(LF) CN lines, respectively. During the observations the spec- 
tral correlator was split in two halves centered at 86.800 GHz 
and 226.767 GHz. This choice allowed us to cover the SiO 
line at 3 mm and the CN lines at 1 mm. The 7=1-0 line 
of H^^CO^ (at 86.754 GHz) is simultaneously covered by the 
3 mm setup. Rest frequencies were corrected for the reces- 
sion velocity initially assumed to be (;o(7/£'L)=1137km s"^ 
The correlator configuration covers a bandwidth of 580 MHz 
for each setup, using four 160 MHz- wide units with an over- 
lap of 20 MHz; this is equivalent to 2000km s"'(770km s"') at 
86.8 GHz (226.8 GHz). Observations were conducted in single 
pointings of sizes 56" and 21" at 3 mm and 1 mm, respectively, 
centered at q'2()oo=02''42'"40.71' and 5200()=-00°00'47.94". The 
latter corresponds to the nominal position of the AGN core, 
as determined from different VLA and VLBI radio continuum 
images of the galaxy (e.g., Gallimore et al. 1 19961 1. Visibilities 
were obtained using on-source integration times of 20 minutes 
framed by short (~2min) phase and amplitude calibrations on 
the nearby quasars 0235+164 and 0336-019. The absolute flux 
scale in our maps was derived to a 10% accuracy based on 
the observations of primary calibrators whose fluxes were de- 
termined from a combined set of measurements obtained at the 
30m telescope and the PdBI aiTay. The bandpass calibration was 
caiTied out using NRAO150 and 1749-1-096 and is accurate to 
better than 5%. 

The image reconstruction was done using standard 
IRAM/GILDAS software (Guilloteau & Lucas 2000). We used 
natural weighting and no taper to generate the SiO line map 
with a size of 64" and 0.25"/pixel sampling; the coiTesponding 
synthesized beam is 3.5" x 2.1", PA=29°. We also used natu- 
ral weighting to generate the CN map with a size of 32" and 
0.12"/pixel sampling; this enables us to achieve a spatial res- 
olution of ~1" (1.5" X 0.8", PA=19°). The conversion factors 
between Jy beam ' and K are 22 K Jy ' beam at 86.8 GHz, and 
21 K Jy~' beam at 226.8 GHz. The point source sensitivities are 
derived from emission-free channels. They are 1.5 mJy beam ' 
in 2.5 MHz-wide channels at 3mm, and 7.6 mJy beam ' in 
6.25 MHz-wide channels at 1mm. Images of the continuum 
emission of the galaxy have been obtained by averaging those 
channels free of line emission at both frequency ranges. The cor- 
responding point source sensitivities for continuum images are 
0.3 mJy beam-' at 86.8 GHz and 0.7 mJy beam^' at 226.8 GHz. 

In this work we use the 1.8" resolution CO(l-O) interfer- 
ometer maps of NGC 1068 obtained with the PdBI by SOO. We 
list in Table 1 the relevant parameters for the SiO and CN PdBI 
observations presented in this paper. 

2.2. Single-dish spectra 

New CH3OH and HNCO observations of the nucleus of 
NGC 1068 were carried out in July 2006 with the IRAM 30m 
telescope at Pico de Veleta (Spain). We observed the 2^-1 ^ group 
of transitions of CH3OH. This group is a blended set of four 
lines. Hereafter we consider the frequency of the 2o-loA^ line 
(96.741 GHz) as velocity reference. We observed the higher fre- 
quency group of transitions of methanol at 2 mm (denoted as 3^- 
2^), a blended set of eight lines. The frequency of the 3()-2oA^ 
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Fig. 1. The continuum maps obtained with the PdBI towards 
the nucleus of NGC 1068 at 3 mm (grey scale contours: 3%, 
5%, 10%, 15%, 25% to 85% in steps of 15% of the maxi- 
mum=39.8 mJy beam"') and 1 mm (white contours: 20% to 
95% in steps of 15% of the maximum=14.2 mJy beam"'). 
The lowest contours coiTespond to ~4cr levels at both frequen- 
cies. (Ac, A5)-offsets are relative to the phase tracking center 
(Q'2000=02''42'"40.71', 520()0=-00°00'47.94"), which coincides 
with the position of the AGN (filled square). Black-filled and 
white-filled ellipses show the beams at 3 mm (3.5"x2.1" at 
PA=29°) and 1 mm (1.5"x0.8" at PA=19°), respectively. 



line (145.103 GHz) is taken as velocity reference. We also ob- 
served the 4o4-3o3 rotational transition of HNCO (hereafter des- 
ignated as 4-3) at 87.925 GHz. The corresponding beam sizes of 
the telescope are ~28" at 3 mm and 17" at 2 mm. The 3 mm and 
2 mm SIS receivers of the 30m telescope were tuned to the red- 
shifted frequencies of the lines around iio(//£'L)=1137km s"'. 
The velocity range covered was 1600 km s"' for the 3 mm lines 
and 900 km s ' for the 2 mm lines. The wobbler switching mode 
was used to obtain flat baselines with a maximum beam throw 
of 4'. 

Typical system temperatures during the observations were 
~200 K at 3 mm and -300 K at 2 mm. All receivers were used 
in single side-band mode (SSB), with a high rejection of the im- 
age band. The calibration accuracy is better than 20%. Pointing 
of the 30m telescope was regularly checked every 1 .5 hours by 
observing a nearby continuum source; we found an average rms 
pointing eiTorof 2-3" during the observations. 

Throughout the paper, line intensities are given in antenna 
temperature scale, T*. The T* scale relates to the main beam 
temperature scale, r^b, by the equation Tmb = (^eff/fieff)r*, 
where Feff and Betf are, respectively, the forward and beam ef- 
ficiencies of the telescope at a given frequency. For the IRAM 
30m telescope Feff/Beff = 1.26 (1.48) at 86 GHz (145 GHz). 

3. Continuum maps 

Figure [T] shows the continuum maps derived at 86.8 GHz and 
226.8 GHz in flie nucleus of NGC 1068. The morphol ogy o f 
the maps is similar to that described by Krips et al. d2006l l. 
who observed the continuum emission at the nearby frequen- 
cies 110 GHz and 230.5 GHz. The emission consists of a central 
component (AGN core), a NE elongation (jet), and a S W elonga- 
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Table 1. Observational parameters of the SiO and CN PdBI data. 





SiO 


CN 


Line 


v=0, J=2-l 


N=2-l (/=5/2-3/2 (HF), 7=3/2-1/2 (LF)) 


Frequency 


86.8 GHz 


226.9 GHz (HF), 226.7 GHz (LF) 


Beam 


3.5"x2A",PA=29° 


1.5"x0.8", PA=19° 


1-cr (rms) 


1.5 mjy beam-i (Av=2.5 MHz) 


7.6 mJy beam"' (Ay=6.25 MHz) 


Field of view 


56" 


21" 


T„,b/S 


22KJy-' beam 


21KJy-' beam 



tion (counter-jet). The counter-jet and the jet are only detected 
at 3 mm. 

We have used the GILDAS task UV-FIT to fit the contin- 
uum visibilities at both frequencies using a set of three point 
sources. Most of the flux comes from the compact source lo- 
cated at the position of the VLBI radio core at (Aa, A(J)~(0", 
0"). The AGA^ core has a flux of 39.8+0.3 mJy and 16.0+0.7 mJy 
at 86.8 GHz and 226.8 GHz, respectively. The NE jef component 
is fitted by a point source of 23.8+0.3 mJy at 86.8 GHz, located 
at (Aa, A5)=(2.6", 3.4"). The 3 mm counter-jet point source at 
(Aa, A5)=(-1.6", -4.0") has a flux of 4.7+0.3 mJy. There are 
nevertheless hints of extended emission, mostly at 3 mm, not 
fully accounted for by the point source fitting (Krips et al. 2006). 
Taking into account the differences in frequency, spatial resolu- 
tion and fitting functions used, the solution described above is 
in agreement with that found by Krips et al. t2006j (see their 
Table 1 for a detailed description). 

4. Molecular gas distribution 

4.1. SiO distribution 

Figure |2] shows the SiO(2-l) intensity map obtained by in- 
tegrating the emission in velocity channels from u = 967 to 
1307 km s"'. In this velocity range line emission is detected 
at significant >3cr levels on scales larger than the beam inside 
the PdBI field of view. We overlay the SiO(2-l) map on the 
CO(l-O) map of SOO. The CO map fairly represents the over- 
all molecular gas distribution in NGC 1068. SOO estimated that 
the CO PdBI map recovers ^80% of the total flux, based on the 
comparison between the CO flux measured by the PdBI and the 
flux derived from the CO maps obtained through the combina- 
tion of the BIMA array and NRAO single-dish data of Heifer & 
Blitz (119951) . For the sake of comparison, the spatial resolution 
of the CO(l-O) map shown in Fig.|2l originally ~1.8", has been 
degraded to that of the SiO data. In Fig.|2]we also compare the 
SiO map to the morphology of the jet derived from the 3 mm 
continuum data. 

Unlike CO, most of the total SiO emission detected inside 
the PdBI primary beam comes from a circum-nuclear molecular 
disk located around the AGN, referred to as CND. The integrated 
flux of SiO in the CND is ~2 Jy kms ' . This is -74% of the to- 
tal flux detected by U04 at 28" spatial resolution with the 30m 
telescope towards the AGN. U04 estimate that ^75% of the SiO 
emission of their AGN spectrum comes from the CND itself, 
with a residual ^25% contribution from the SB ring. Taken to- 
gether, these estimates indicate that the PdBI recovers practically 
all of the SiO flux in the CND. Little SiO emission is detected 
from the SB ring in the PdBI map. A few SiO clumps, show- 



ing sizes ~beam and intensities at the 3-4-cr level, are distributed 
along the SB ring. 

The overall distribution of SiO in the CND, spatially re- 
solved by the PdBI beam, shows an east-west elongation. The 
deconvolved full size of the emission at a 4cr level is ~400 pc 
along the east-west axis. The CND is marginally resolved along 
the north-south axis; the corresponding deconvolved size is 
<150 pc. These sizes are comparable with those derived from the 
higher resolution 1-0 and 2-1 CO maps of the CND published 
by SOO, as well as with those derived from the CN(2-1) map 
described in Sect. 14.21 The SiO emission is slightly asymmetric 
around the AGN: the SiO emission peak is at (Aa, A5)~(0.5", 
0"). Instead, the CO map shows a maximum at (Aa, A6)~{1.2", 
0") (see Fig. 1 of SOO). The latter corresponds to the E CO knot, 
in the notation of SOO. The E CO knot is thus shifted ~50 pc east 
relative to the E SiO knot. 



4.2. CN distribution 

Both lines of CN (LF and HF, as defined in Sect. |2]i are detected 
in the CND of NGC 1068. As expected, none of the CN lines 
are detected towards the SB ring, as this region is located well 
beyond the edge of the primary beam of the PdBI at 226.8 GHz. 
The HF/LF intensity ratio shows no significant spatial variation 
inside the CND. The average ratio, ~1.7, is in close agreement 
with the value expected in the limit of optically thin emission 
(HF/LF~ 1.64 + 0.14). In the remaining of this paper, we will 
use the higher S/N ratio maps of the HF line. 

Figure [3] shows the intensity map of the HF line of CN(2- 
1) obtained towards the CND of NGC 1068. We have integrated 
the emission in velocity channels from y = 955 to 1277 km s"', 
i.e., similar to the velocity range used to derive the SiO map. 
With this choice we encompass all significant line emission and 
at the same time avoid blending with the LF CN(2-1) line. At 
this spatial resolution {~l"), the overall morphology of the CND 
traced by CN roughly resembles that seen in CO lines. At close 
sight there are noticeable differences between these tracers, how- 
ever. The two CN knots, connected by a bridge of emission on 
the northern side of the CND, lie at (Aa, A6)~{0.9'\ 0") and 
(Aa, A6)~(-l.7", -0.5"). As for SiO, the CN knot located east 
is significantly closer to the AGN compared to the corresponding 
CO knot. 

5. Molecular gas kinematics 

5. 1 . Background from previous work 

The overall kinematics of molecular gas in the disk of NGC 1 068 
were modeled by SOO and BOO. One of the scenarios advanced 
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Fig. 2. a) (Upper panel) The SiO integrated inten- 
sity map (contour levels are 3cr to 12cr in steps of 
lo-=0.082 Jy km s"' beam"') is overlaid on the CO(l-O) 
integrated intensity map of SOO (color scale as shown in units 
of Jy km s"' beam"'), b) (Middle panel) The same as a) but 
showing a zoomed view on the inner 12" around the AGN 
(identified by the cross), c) (Lower panel) The SiO map con- 
tours are overlaid on the 3 mm continuum map (color scale) of 
Fig.[T] The SiO beam is shown in aU the panels by an ellipse. 




Aa(arcsec) 



Fig. 3. The CN integrated intensity map. Contour levels 
are 3cr, 4cr, 6cr, 9cr, 13cr to 25cr in steps of 6cr, with 
lcr=0.40 Jy km s"' beam"'. The AGN is identified by the cross 
and the CN beam is shown by an ellipse. 



by SOO invokes the presence of two embedded bars: an outer 
oval and a nuclear bar of ~17 kpc and 2.5 kpc deprojected diam- 
eters, respectively. The two bars are coupled through resonance 
overlapping, so that the corotation of the nuclear bar coincides 
with the outer inner Lindblad resonance (ILR) of the oval. The 
gas response inside r ~5"(350 pc), a region that encompasses 
the whole CND detected in SiO and CN, would lie between the 
ILRs of the nuclear bar. 

SOO and BOO proposed also a nuclear warp model to alterna- 
tively account for the anomalous kinematics of molecular gas in 
the CND. The non-coplanar gas response in the nuclear region 
could have been caused by the interaction of the molecular disk 
with the jet or with the associated ionization cone. BOO's model 
proposes a warped molecular disk that evolves from a low incli- 
nation (/ -30°) at r -2.5" to a highly-inclined disk (/ -77°) at 
r -0.5". This configuration naturally explains the steep velocity 
gradient (high velocities) measured in the major axis position- 
velocity (p-v) diagram of CO at small radii (r <2"). In particu- 
lar, SOO favor a hybrid model where the gas response to the bar is 
combined with a nuclear warp. Besides being able to reasonably 
fit the gas kinematics, the main advantage of the hybrid solution 
over the pure bar model is that the first explains the asymmet- 
ric excitation of molecular gas evidenced by the remarkably dif- 
ferent R-2-1/1-Q line ratios measured in the east and the west 
CO lobes. BOO argues that the east lobe, the one showing the 
higher R, corresponds to the directly X-ray illuminated surface 
of a warped disk, whereas the west lobe corresponds to the back 
of the warped disk, characterized by a lower R. The illumina- 
tion of the CND gas by X-rays may also explain the particular 
chemistry of molecular clouds analyzed in Sect. 16.21 

The kinematics of molecular gas in the inner r -50 pc 
region are also complex, as recently revealed by the 2.12//m 
H2 1-0 S(l) map of Muller-Sanchez et al. (|2009] |. These data, 
mostly sensitive to hot (T^^ ^10^ K) and moderately dense 
(n(H2)-10-' cm""*) molecular gas, show the existence of a struc- 
ture bridging the CND and the central engine. This connection 
is made through highly elliptical streamers detected in H2 lines. 
Davies et al. (,2008) used these data to study the large-scale kine- 
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matics of molecular gas in the CND itself and concluded that 
the whole CND is a lopsided ring in expansion (Krips et al., 
in prep.). The scenario of an expanding ring, first advanced by 
Galliano & AUoin (I2002I I. could be Unked to the onset of a nu- 
clear warp instability. 



5.2. SiO kinematics 

Figure HJz shows the isovelocity contour map of the CND of 
NGC 1068 derived from SiO. Isovelocities have been obtained 
using a 3cr clipping on the data cube. The velocity centroid of 
SiO emission towards the AGN is uacn ~1122+10 km s"' in 
HEL scale. Gas velocities are red {u > uacn) on the northwest 
side of the CND, while blue velocities (u < vacn) appear south- 
east. This picture is consistent with a spatially resolved rotat- 
ing molecular structure. The kinematic major axis is fitted to 
PA ~105+10° to expand the maximum range of SiO radial ve- 
locities at the edge of the CND (r ^2.5")- This solution lies 
within the wide range of PA values derived from the photomet- 
ric fitting of the disk and, also, from the observed stellar and 
ga s kine matics (PA ~80°-125°; Br inks e t al. [19971 Dehnen et 
al. [19971 SOO; BOO; Emsellem et al. [20061 Gerssen et al. 2006). 
We note that SiO isovelocities are twisted to larger position an- 
gles {PA ~135°) at smaller radii (r <1"), which already suggests 
the presence of non-circular and/or non-coplanar motions in the 
CND. A similar isovelocity twist is seen in other molecular trac- 
ers, as pointed out by BOO. Outside the CND, the large-scale 
stellar and gas kinematics indicate that the major axis of the disk 
is oriented east-west (Emsellem et al. .2006j . We therefore adopt 
in the following PA=90° as the best guess for the major axis ori- 
entation. 

Figure Hf? shows the SiO p-v diagram along the major axis. 
The p-v diagram is mostly symmetric for the lowest contour lev- 
els around the position of the AGN (Ax' =0"), and with respect 
to a radial velocity t;=i;5^.5(HEL)=1140+10 km s"'. We there- 
fore adopt the first as the dynamical center and the second as 
the systemic velocity Vsys of the galaxy. This value of Vsys is in 
good agreement with all previous determinations of the galaxy's 
receding velocity (~1137+3 km s"'; e.g., de Vaucouleurs et 
al.[T991 ; Huchra et al. 1999). The highest contour levels of the p- 
V plot are shifted east, however; this reflects the overall east-west 
asymmetry of the CND described in Sect l4.1l 

The steep velocity gradient {high velocities) measured along 
the major axis at r <2" can be interpreted at face value as a 
signature of rotation for a highly inclined disk. Fig.|4j7 neverthe- 
less confirms the presence of non-circular motions in the CND: 
SiO emission is detected at regions of the p-v diagram which are 
forbidden by circular rotation: emission in quadrants I and II of 
Figure llf?, can be interpreted by gas being in apparent counter- 
rotation and/or moving radially outward. In either case, this can 
be taken as evidence that gas orbits must be either elongated, 
if they lie in a common plane, or, alternatively, non-coplanar Q. 
Note that the detection of SiO gas in these regions cannot be ex- 
plained by beam smearing effects: the deconvolved size of the 
emission for the velocity channels showing the largest excur- 
sions into I and II (~1155 km s"' and 1 125 km s"', respectively) 
is ^1.5x the beam size at PA ~90°. Similar features are detected 
in the 0.7" resolution CO(2-l) map of SOO (Fig.4 of SOO). As 
discussed in Sect. 15.31 CN emission is also associated with/or- 
bidden velocities. 



' in this context a combination of elongated and non-coplanar orbits 
cannot be excluded. 



In the extended ILR region of the CND it is expected that 
the gas flow would trace the gradual transition from the xi orbits 
of the nuclear bar (between corotation and the outer ILR) to the 
X2 orbits (between the ILRs) (e.g., Athanassoula 1 1 9921 Buta & 
Combes , 1996J . Orbits of the xi family would thus account for 
most of the CO and SiO emission detected at high velocities at 
r <2" in the major axis p-v diagrams of Fig.|5]3. Contrary to SiO 
emission, which is mainly restricted to the CND, CO emission 
is detected over the SB ring (Ax' =+14") but also in a region 
that connects the CND with the SB ring. A significant percent- 
age of the CO emission on these intermediate scales (Ax'=+3- 
8") is detected at low velocities, i.e., at velocities significantly 
below the terminal limit imposed by the rotation curve at this 
radii (Fig. [Sj?). This is the kinematic signature of the gas re- 
sponse to the stellar potential, characterized in this region by the 
crowding of x\ orbits at the leading edges of the bar (SOO; BOO). 
As expected, the x\ + xi orbit combination produces a double- 
peaked line-of-sight velocity distribution in the CO p-v diagram 
that shows the characteristic figure-of-eight trend as a function 
of radius where high velocities lie at small radii (Kuijken & 
Merrifield [T995l Garcia-BuriUo & GueUn 1995). 

In the light of the published dynamical models of NGC 1068 
(SOO; BOO), we conclude that most of the SiO emission in the 
CND is associated with a region identified as the extended ILR 
region of the galaxy. Non-coplanar motions have alternatively 
been proposed to fit the CND kinematics. Nonetheless in either 
case (nuclear bar or nuclear warp) the high velocity emission 
detected in SiO is physically associated with gas lying at small 
radii (r <2"). 



5.3. CN kinematics 

Figure |6]3 shows the isovelocity contour map of the CND of 
NGC 1068 derived from CN. Isovelocities have been obtained 
using a 3cr clipping on the data. Compared to SiO, the higher 
spatial resolution of the CN maps provides a sharper but more 
complex picture of the gas dynamics in the CND. The rotating 
pattern of the CN emitting gas is highly perturbed. The overall 
kinematics suggest an east-west orientation of the major axis. 
However, at close sight the PA of the major axis changes from 
90° to 180° over the spatial extent of the E CN knot. In addi- 
tion, isovelocities show an irregular pattern at the W CN knot. In 
spite of the more complex kinematics revealed by the higher res- 
olution CN map on smaller scales, we will nevertheless assume 
in the following the same overall orientation and kinematic pa- 
rameters determined from SiO. As argued in Sect. 15.21 this is a 
reasonable guess compatible with previous determinations based 
on the observed stellar and gas kinematics. 

Figure |6l7 shows the CN p-v diagram along PA ~90°. At 
this spatial resolution the p-v diagram shows a strong east-west 
asymmetry: most of the gas emission appears associated to the E 
knot. Similarly to SiO, CN emission is detected at regions of the 
p-v diagram which are forbidden by circular rotation. For CN, 
forbidden velocities appear mostly at quadrant I (see Fig. jSJ?). 
Due to the higher spatial resolution, the velocity gradient mea- 
sured along the major axis is ~a factor of two steeper in CN than 
in SiO. Otherwise, the two molecular tracers of the CND share 
the same kinematic features, such as the existence of high ve- 
locity emission at small radii (for CN r <!"; see Fig.|6j7), likely 
connected to X2 orbits. Like SiO, CN shows no detectable emis- 
sion at low velocities on intermediate scales (r >3"; see Fig.|5]/7). 
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Fig. 4. a) (Upper panel) SiO isovelocities contoured over false- 
color velocity maps. Velocities span the range [1050 km s"', 
1225 km s"'] in steps of 25 km s"^ Velocity scale is HEL. The 
AGN position is marked with a star. The position of the kine- 
matic major axis derived from SiO at PA ~105° is indicated by 
the dashed line. An ellipse shows the beam, b) [Lower panel) 
Position-velocity diagram of SiO along PA=90°. Contour lev- 
els are -3, 3, 3.75, 4.5 to 9 in steps of 1.5 mJy beam"'. The 
position of the AGN (Ax'=0") and Vsys(^EL)-ll40 km s"' are 
highlighted. We indicate the position of quadrants I and II dis- 
cussed in the text. 
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Fig. 5. a) (Upper panel) Position-velocity diagram of CO(l-O) 
(contours; data from SOO degraded to the resolution of SiO) and 
SiO(2-l) (color scale; this work) along PA=90°. We highlight 
the kinematic signature of the bar leading edges, formed by the 
precession of xi orbits, seen in the low velocity CO component. 
We plot the projected terminal velocities (empty squares) de- 
rived from the rotation curve fitted by SOO. b) (Lower panel) 
Same as a) but showing the position-velocity diagram of CO(l- 
0) (contours; data from SOO) and CN(2-1) (color scale; this 
work), degraded to the resolution of CO, along PA=90°. 



6. Molecular abundances 



In Sect. Owe use the average SiO/CO and SiO/H'3CO+ inten- 
sity ratios measured in the CND and the SB ring to derive the 
global abundance of SiO in the two regions, using a one-phase 
LVG radiative transfer model. We further analyze in Sect. l6.2l the 
SiO/CO and CN/CO intensity ratio patterns inside the CND, as 
derived from the high-resolution PdBI maps. We discuss how 
these patterns can be interpreted in terms of different physi- 
cal/chemical properties of molecular gas in the CND. 



6.1. SiO abundances in the CND and the SB ring 

We have estimated the average SiO/CO velocity-integrated in- 
tensity ratios (< /?sio/co >) measured in the CND and the SB 
ring regions defined as follows. Based on the CO distribution 
shown in Fig.|2]3, we delimit the SB ring as the region between 
r ~10" and r ~20" with SiO intensities >2.5cr levels to derive 
integrated values. Similarly, the CND region is defined using a 
2.5cr clipping on the integrated intensities of the SiO map of 
Fig. I2J7. Intensity units used to derive line ratios are K km s"' 
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LVG results for SiO at Tk = 50 K 
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Fig. 6. Same as Fig. |4] showing in a) (Upper panel) CN isove- 
locities, derived from the high frequency (HF) 2-1 transition of 
CN, contoured over false-color velocity maps. Velocities span 
the range [1025 km s"', 1200 km s"'] in steps of 25 km s"'. 
Similarly, we show in b) (Lower panel) the p-v diagram of CN 
along the major axis at PA=90°. Contour levels are -15.2, 15.2, 
22.8, to 68.4 in steps of 7.6 mJy beam"'. Velocities are referred 
to the HF 2-1 line. We highlight the range of the HF and LF lines 
ofCN. 



with K in T^b scale. The < /?sio/co > ratios are thus equivalent 
to brightness temperature ratios, as the line widths for SiO and 
CO are comparable at this spatial resolution. 

The maximum ratio in the SB ring, corrected by the different 
primary beam attenuation factors, is ~0. 004+0. 001, a ratio esti- 
mated from the average spectrum of the clumps detected in SiO. 
This is a factor of ~20 lower than the corresponding ratio av- 
eraged over the CND, estimated as ~0. 08+0. 01. The dichotomy 
between the CND and the SB ring is apparent from the remark- 
ably different < /?sio/co > ratios measured in the two regions. 
The origin of such a dichotomy could be found in the differ- 
ent physical conditions or, alternatively, in the different chemical 
properties of molecular gas in these two regions. In the follow- 
ing, we use a radiative transfer code to fit the set of available line 
ratios in order to disentangle the contribution of different factors. 
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Fig. 7. LVG results for SiO, assuming a kinetic temperature of 
Jk = 50 K. The x and y axes of the panel are, respectively, 
the SiO column density per velocity width (A^sio/Ay) and the 
volume density of H2 molecules (n(H2)) in a cloud. The solid and 
dashed lines indicate constant SiO(3-2)/SiO(2-l)s< /?32/2i > 
line ratio and SiO(2-l) brightness temperature, respectively. The 
tags "SB" and "CND" identify the mean values measured in the 
SB ring and the CND. The < R32/21 > ratios within Icr from 
the mean values in the CND and the SB correspond to the light- 
and medium-gray strips. The two strips overlap in a thinner band 
painted in darker gray. 



6.1.1. LVG model results 

To estimate the SiO abundances both in the CND and in the SB 
ring, we adopt a one-phase LVG approach to model the radia- 
tive transfer of the SiO emission. We use an LVG code to trans- 
late the SiO brightness temperatures into the SiO column density 
per velocity width (Nsio/Au) and the volume density of its colli- 
sional partner, H2 (n(H2)), in each cloud. We feed into the LVG 
code the average SiO(2-l) brightness temperature (from our 
PdBI observations) and the average SiO(3-2)/SiO(2-l) velocity- 
integrated line ratio (from U04), hereafter < 7^32/21 >, measured 
in each target region. Provided that the 3-2 and 2-1 SiO emis- 
sions arise from the same clouds in each region, < 7^32/21 > re- 
mains independent of the unknown filling factor rjf. To further 
constrain the model, the gas kinetic temperature is fixed at 50 K, 
as obtained by Tacconi et al. ( 1994) from multi-transition CO ob- 
servations of the CND. This assumption is not critical, since the 
results barely depend on the gas temperature below, e.g., 100 K. 
Quite likely, 100 K is a comfortable upper limit to the average 
temperature of the SiO emitting gas in a PdBI beam. 

The output of the LVG code is summarized in Fig.|7] Since 
the filling factors are unknown the possible solutions for the 
CND (SB ring) clouds would lie over the corresponding < 
^32/21 > isocontour, anywhere to the right of the intersect with 
the CND (SB ring) SiO(2-l) isocontour (i.e., the lower rjf, the 
higher A^sio/Au in a cloud). The loci of possible solutions be- 
come two-dimensional bands, once we also consider the uncer- 
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Fig. 8. SiO abundances estimated from the LVG model as a func- 
tion of the SiO column density per velocity width (Nsio/Au) in a 
cloud (see text for details). Each blue, black or red curve is de- 
fined by three parameters: < /?32/2i >, < ^sio/co > and X. The 
conditions in the SB ring are represented by the red curves and 
those in the CND by the others. The solid curves correspond to 
the measured < 7^32/21 > and the dotted ones to this value +0- 
level. The SB solutions are only given for a Galactic X, while 
in the case of the CND a five times lower X is also considered 
(blue). The areas where the most likely solutions lie are hatched 
as explained in the text. They range in Nsio/Au they cover is so 
that the CO(l-O) brightness temperature in a cloud is neither 
smaller than the measured brightness temperature, nor higher 
than the gas kinetic temperature (50 K). 



spond to typical gas densities of ~ 10^ cm ^ in the CND, while 
they are > 5 times lower in the SB ring. 

We search for a reasonable estimate of the SiO abundance 
with respect to H2 (X(SiO)= Nsio/Nn.) in the SB ring and the 
CND. For this we use ancillary CO observations of NGC 1068 
to infer the H2 column density. A full radiative transfer ap- 
proach (e.g., LVG) is not possible, however, in the SB ring, 
as multi-transition CO observations are not available in this 
region. Alternatively, we use the standard conversion between 
the CO(l-O) integrated intensities and the H2 column densities 
(A^H,) given by the so called X factor in the two target regions. 
The X conversion factor is of order ~ 2 x 10^° cm"^/(km s"') 
in the Milky Way (Strong et al. 1988), a value which roughly 
holds true for most nearby galaxies. For the SB ring, we adopt 
the Galactic X value. It is known, however, that a lower X has 
been found in galaxy centers, where molecular clouds are sub- 
ject to abnormal conditions (e.g., they are unlikely to be virial- 
ized). In particular, X in the CND of NGC 1068, could be ~ 5 
times lower than the standard value, according to the estimates 
of U04, a result which is in agreement with the findings obtained 
in the nuclear disks of other spiral galaxies (e.g., Garcia-Burillo 
et al. 1993). 

Under these assumptions, we proceed as follows to estimate 
X(SiO) in the SB ring and the CND: 



- Each target region is characterized by fixed < 7^32/21 > and < 
^sio/co > values Q assumed to be independent of the filling 
factor, and by a certain X. 

- The measured < /?32/2i > allows us to calculate the SiO(2- 
1) brightness temperature as a function of A^sio/Af from the 
diagram in Fig.|7l 

- Ah, /Ay is proportional to the product of < .Rsio/co > and the 
SiO(2-l) brightness temperature. Thus, ultimately, X(SiO) 
can be expressed in terms of a certain function of Nsio/Av 
parameterized by < /?32/2i >, < ^sio/co > ratios and X. 

- The range of values for Nsio/Au can be constrained. The 
lower limit to Asio/Au corresponds to rjf - 1. On the other 
hand, for rjf arbitrarily low, the CO(l-O) brightness temper- 
atures in a cloud would rise above the assumed kinetic tem- 
perature (50 K), which cannot occur if CO is coUisionally 
excited. This sets an upper limit to Asio/Ay. 

The end product of this analysis is the plot in Fig. [8] There 
we represent X(SiO) as a function of Nsio/Au for the clouds in 
the SB ring and the CND. In the CND we explore two solu- 
tions with X equal to the standard factor and to a five times 
lower value. In either case, we note that the values of X(SiO) 
in the CND, ~(l-5)xlO"^, are significantly boosted with respect 
to those typically found in Galactic dark clouds (< 10 '^; e.g., 
Ziurys et al.'1989); this result calls for a mechanism of efficient 
dust grain processing in the CND. It is worth noting that the true 
SiO abundances could be even higher than our estimates if only a 
fraction of the CO emitting clouds were simultaneously emitting 
in SiO (the curves in Fig. [8] would have to be shifted upwards). 
We also note that X(SiO) is at least >one order of magnitude 
higher in the CND than in the SB ring. Actually, were X ~5 
times lower in the CND than in the SB ring, the SiO abundances 
in the two target regions would differ by up to ~two orders of 
magnitude. 



6.1.2. SiO/H'3CO+ ratios 

The results of the radiative transfer model described above can 
be compared to an alternative estimate of the abundance of SiO 
in the CND, based on the observed SiO/H'^CO+ fine ratio. The 
H'^CO^(l-O) line, observed simultaneously with the SiO(2-l) 
line, has also been detected in the CND of NGC 1068. Compared 
to SiO, the H'^CO^ line is ~a factor of three weaker. This ratio 
is significantly larger than that derived by U04 who reported a 
value ~1 from a 30m spectrum of the CND. The disagreement 
between the PdBI and the 30m results is an indication that, com- 
pared to SiO, the contamination from the SB ring in H'^CO^ is 
more severe, and cannot be accurately estimated using the lower 
resolution 30m map. 

Given the weakness of H'^CO^ in the CND, we omit any dis- 
cussion on the 2D morphology of the emission of this line. We 
have nevertheless used H'^CO^ to independently estimate the 
abundance of SiO in the CND using a LVG code to fit the global 
SiO(2-l)/H'-'CO^(l-0) intensity ratio measured in the nucleus 
of NGC 1068. The value of X(SiO) can be derived assuming a 
plausible range of physical conditions for the gas. This relies on 
two basic assumptions (see also discussion in Garcia-Burillo et 
al.EOOOiand Usero et al. l2006] l. First, since SiO and Hi3CO+ 
have similar dipole moments and the observed transitions have 
comparable upper state energies, it is plausible to assume similar 
physical conditions for the two species. In our case we take the 
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Fig. 9. a) (Upper left panel) The SiO(2-l)/CO(l-0) velocity-integrated intensity ratio (/?sio/co) (color scale) derived at the spatial 
resolution of SiO inside the CND. This ratio is equivalent to a brightness temperature ratio for the observed common line widths 
of SiO and CO. An ellipse shows the SiO beam, b) {Lower left panel) The SiO/CO ratio is here displayed as function of velocity 
and spatial offset across the major axis at PA-9Qi° (SiO/CO pv). The CO pv plot is displayed in contours (CO pv). c) {Upper right 
panel) Same as a) but showing here the CN(2-1)/CO(1-0) velocity-integrated intensity ratio (/?cN/co)(color scale) derived at the 
spatial resolution of CO. An ellipse shows the CO beam, d) {Lower right panel) Same as b) but adapted to compare here the CN/CO 
ratio along the major axis (CN/CO pv) with the corresponding CO pv plot. 



CND values estimated by U04 from the multi-line 30m survey 
of NGC 1068: n(H2)~10^^cm-^ and Tk = 50 K. As a second as- 
sumption, we take a standard abundance for H'^CO^=2.5xlO "* 
and derive X(SiO) from the SiO/H'^CO^ column density ratio 
fitted by LVG. In our case this corresponds to a typical abun- 
dance of the main isotope X(H'^CO^) =10 **, and to an isotopic 
ratio ['^C]/['^C]- 40, similar to the values found in the CND 
of our Galaxy (Wannier |1980l l. The derived abundance of SiO 
in the CND is ~2.2xl0"^. This estimate is comfortingly close 
to the lower limit of the X(SiO) value derived from the LVG 
modeling of the SiO/CO ratios described above. This agreement 
confirms that the one-phase LVG approach adopted in Sect. l6.1.TI 



provides a reasonable estimate of the SiO molecular abundances 
in NGC 1068, in spite of the inherent limitations of such a sim- 
plified description of the molecular ISM in this model. 



We thus conclude that the dichotomy between the CND and 
the SB ring, regarding the measured < /?sio/co > and < R-nin > 
ratios, cannot be attributed to a comparatively higher volume 
density of SiO gas in the CND, but mainly to a significant en- 
hancement of the SiO abundance in this region. 
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Fig. 10. a) (Upper panel) SiO abundance as a function of gas 
density in the high velocity gas at the position of the AGN (soHd 
line) and in the low velocity at the edge of the CND (dashed 
Une). The abundances are derived from the mean SiO and CO 
brightness temperatures within the respective velocity windows 
as defined in text. The SiO abundance in the two regions can be 
the same only within the grey area, b) (Lower panel) Same as a) 
but particularized for the density and abundance ranges for CN. 



6.2. Chemical differentiation in the CND 

6.2.1. SiO abundances 

Figure |9]3 shows the SiO(2-l)/CO(l-0) velocity-integrated in- 
tensity ratio (Rsio/co) in the CND, derived using a spatial reso- 
lution corresponding to that of SiO. The changes in /?sio/co are 
more significant along the major axis, where the spatial resolu- 
tion is the highest. /?sio/co changes by ~a factor of three inside 
the CND: it goes from a peak value of ~0. 11+0. 01 in a region 
close to the AGN to a minimum of ~0. 04+0. 02 at the eastern 
boundary of the CND (Aa — 1-3"); ^sio/co reaches a secondary 
maximum (~0. 10+0.01) at the western edge of the CND (Aa — 
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Fig. 11. Spectra obtained with the IRAM 30m telescope in se- 
lected molecular line tracers of shock chemistry. We show in 
a) (Upper panel) the spectra of the 2-1 line of CH3OH. Intensity 
scale is in T* (K) units; velocities are in HEL scale (km s"'). 
The Gaussian fit to the spectrum is displayed in blue, b) (Middle 
panel) Same as a) but for the 3-2 line of CH3OH. c) (Lower 
panel) Same as a) but for the 4-3 line of HNCO. 



3"). This range of values confirms and significantly expands the 
differences in the SiO/CO ratios between the east (~0.07) and the 
west (~0.10) lobes of the CND, identified by U04 as coming, re- 
spectively, from the blue and the red velocity components of the 
line profiles. Figure. |9j7 shows the SiO/CO ratio displayed here 
as a function of velocity and position across the major axis of 
the galaxy. The differences in the SiO/CO ratios, which go from 
~0.05 to ~0.5, expand out to ~one order of magnitude depend- 
ing on the velocity channel and position along the major axis. 
Line ratios follow a regular pattern, where the highest ratios are 
associated with high velocities. In particular, the SiO/CO ratio is 
-0.35+0.05 at I y - y,^, |=80-160km s"' towards the AGN and 
the western secondary maximum. By contrast, the correspond- 
ing average ratios at low velocities (\ v - Usys |<+80 km s"') are 
~a factor of 6-7 lower (-0.05+0.01) at the edges of the CND 
(Ajc'=+3"). 
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The dichotomy between the SiO/CO ratios measured at high 
and low velocities can be interpreted in terms of different phys- 
ical conditions (i.e., different n(H2)), or, alternatively, as a sig- 
nature of chemical differentiation (i.e., different X(SiO)) in the 
CND. Figure. [TOb shows the range of possible LVG solutions 
that fit the CND ratios. We explore two extreme solutions: in 
the first case, a common X(SiO) is assumed, whereas in the 
second case, we adopt a common n(H2). While physical con- 
ditions and chemical abundances in the gas are to some extent 
interconnected, by exploring the two extreme cases described 
above, we can investigate which is the main driving cause of 
the reported changes in the line ratios. We restrict H2 densities 
to lie within a range of values around the average density de- 
rived from the fit of the < /?32/2i > ratio measured in the CND 
(n(H2)=10^ '^''^cm"-'). With these restrictions. Fig. [TOb shows 
that we need >one order of magnitude increase in n(H2) to fit 
the progression of /?sio/co values observed from low velocities 
(at larger radii) to high velocities (at smaller radii). Alternatively, 
the scenario of a common n(H2) but different X(SiO), suggests a 
much wider range of solutions for the CND. In this case, X(SiO) 
would need to be ~a factor of 5-6 higher at high velocities. 
Although a tenfold increase of n(H2) cannot be formally ex- 
cluded based on these datqj, we note that this solution implies 
that the change should occur on scales of ~200-300 pc (~the 
CND radius). This is at odds with the significantly smaller differ- 
ences (~a factor of 2-5) that exist between the average densities 
of SiO gas detected in the SB ring and the CND, two regions 
more than 1 kpc apart. Therefore we can interpret the spatial 
variations of < /?sio/co > in the CND as the likely signature of 
chemical differentiation. In this scenario, the abundance of SiO 
is significantly enhanced at high velocities (i.e., at small radii 
r <2"; see Sect. 15.21) . This interpretation is supported by the cor- 
relation found between < /?sio/co > and the X-ray irradiation of 
the CND discussed in Sect.lT21 

6.2.2. CN abundances 

Figure |9j; shows the CN(2-1)/CO(1-0) velocity-integrated in- 
tensity ratio (^cN/co) derived in the CND, using a spatial reso- 
lution corresponding to that of the CO map. The average ^cn/co 
ratio in the CND is ~0.15. The changes in /?cn/co shown in 
Fig.|9]: are significant along the major axis. /?cn/co changes by 
~a factor of three inside the CND. Similarly to the SiO/CO ra- 
tio, the maximum CN/CO value (~0.3±0.03) corresponds to a 
region very close to the AGN. ^cn/co decreases to <0. 1+0.04 
at the edge of the CND. Figure. |9li shows the CN/CO ratio here 
displayed as a function of velocity and position across the ma- 
jor axis of the galaxy. Similarly to the case of SiO discussed 
above, an inspection of Fig. |9j/ shows that the differences in 
the measured CN/CO ratios are boosted to reach ~one order 
of magnitude depending on the velocity channel and position 
along the major axis. Line ratios follow a regular pattern, with 
the highest ratios being associated with high velocities. Towards 
the AGN, the highest CN/CO values (-0.7+0.09) correspond to 
the high velocity component of the emission. This is in clear 
contrast with the corresponding values measured at low veloci- 
ties (~0. 1+0.02). Following the scheme of Sect. I6.2.T1 we have 
explored two families of LVG solutions that fit the range of 
CN/CO ratios of the CND assuming a common X(CN) or a 
common n(H2). For simplicity we restrict H2 densities to lie 
within the same range values explored for SiO. Perez-Beaupuits 



et al. (120091 1. based on a multi-transition CN analysis, have de- 
rived a more restricted density range (n(H2)=10^'^'cm"^), which 
nevertheless lies within the range explored here. In the first sce- 
nario we need about an order of magnitude increase in n(H2) to 
fit the progression of ratios from low-velocities to high-velocities 
(see Fig.fTOb). The corresponding range of solutions in the sec- 
ond scenario requires ~a factor of 5-6 increase in X(CN). For 
reasons similar to those described in Sect. 16.2.11 we favor an in- 
terpretation of the spatial variations of < ^cn/co > in the CND in 
terms of chemical differentiation. The abundance of CN would 
be significantly enhanced at high velocities (i.e., at small radii 
r <1"; see Sect. 15.31 1. This interpretation is also supported by the 
correlation found between < ^cn/co > and the X-ray irradiation 
of the CND discussed in Sect.lTTil 



' A high-resolution < ^^32/21 > map would be required to quantita- 
tively probe SiO density changes inside the CND. 



7. Molecular gas chemistry in the CND of NGC 1068 

In the following we consider the pros and cons of two differ- 
ent scenarios, shocks and XDR chemistry, regarding their abil- 
ity to explain the molecular abundances measured in the CND 
of NGC 1068 for SiO and CN. Complementary information 
provided by other molecular tracers is a key for this discus- 
sion. Furthermore, we use the information extracted from the 
observed kinematics of the gas and the relation of the derived 
molecular abundances with the X ray irradiation in the CND. 



7. 1 . Shock chemistry 

To explore the prevalence of shock chemistry in the CND of 
NGC 1068 we have used the IRAM 30m telescope to observe 
a set of lines of two molecular tracers of shocks: CH3OH and 
HNCO. We discuss below the line ratios derived for CH3OH, 
HNCO, SiO and CN in NGC 1068, and compare these with the 
ratios derived in galactic and extragalactic templates of shock 
chemistry. We also analyze the potential drivers of shocks in the 
CND. 



7.1 .1 . Tracers of shocks: CH3OH and HNCO 

On theoretical grounds, large abundances of methanol can only 
be produced in gas-phase via evaporation and/or disruption of 
icy mantles. Fast shocks (vsiwck >15-20 km s"') can destroy 
the grain cores, liberating refractory elements, like Si, to the 
gas phase (Caselli et al. [19971 Schilke et al. [T997b . By con- 
trast, slow shocks (vshock <10-15 km s'^) are able to process 
the icy grain mantles, but not the grain cores (Millar et al. |1991t 
Charnley et al. 1995). The different location of Si-bearing ma- 
terial (cores) and solid-phase CH3OH (mantles) in dust grains 
implies that SiO and CH3OH are good tracers of fast and slow 
shocks, respectively. It is expected that the disruption of dust 
grains by slow shocks can inject grain mantle material into 
the molecular ISM without destroying the molecules (Bergin 
et al. 11998! Martin-Pintado et al .200rj . In particular, the abun- 
dance of methanol is seen to be enhanced in shocks associated 
with molecular outflows by more than two orders of magnitude 
over the values typically derived in cold molecular clouds (e.g., 
Bachiller & Perez-Gutierrez 1997; Perez-Gutierrez 1999). In ex- 
ternal galaxies Meier & Turner (2005) associate over-luminous 
methanol lines to shocks in IC 342. The close association be- 
tween SiO and CH3OH emission in IC 342, discussed by Usero 
et al. (2006) corroborates this picture. 

Although it is still debated what is the main production 
mechanism of HNCO, models involving dust grain chemistry 
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Table 2. Line ratios measured in prototypical regions. 



Source 


CH, 0/7(2-1) 
SiO(2-l) 


CHj,0H(i-2) 
SiO(3-2) 


HNCO(4-3) 
SiO(2-l) 


5iO(2-l) 
CN(2-\) 


CN(2-1) 
HCN( 1-0) 


References 


NGC1068-CND 


2.3* (1.8) 


2.1* (1.7) 


1.2* (1.2) 


0.9 


0.25 


1,2 


L1157-B1 
L1157-B2 


2* (1.7) 
2.3* (2.1) 


1.9* (1.7) 
2.9* (2.7) 


- 


3* (7.5) 
>7* (>12) 


0.11* (0.04) 
<0.04* (<0.03) 


3,4 
3,4 


IC342-Ndisk 
IC342-N arm 


5 


4.5* (3.6) 


3.5 


- 


- 


5,6 

7,8 



Notes. We list line ratios derived in the CND of NGC 1068 (NGC1068-CND), the Bl and Bl bullets of the young stellar object LI 157, as well as 
the nuclear disk and Northern spiral arm of the nearby galaxy IC342 (IC342-Ndisk and IC342-N arm, respectively). Line ratios are given in T,,,/, 
units; these are corrected by source coupling factors for single-dish observations. Uncorrected ratios are identified by an asterisk; corrected ratios 
(in Tb-source units) are given in parentheses when necessary. 



References. (1) U04; (2) this work; (3) N isini e t al. i2007l ; (4) Perez- Gutierrez (1999) ; (5) Huttemeister et al. fT997] l; (6) Garci'a-Burillo & 
Martin-Pintado ( l200Tt ; (7) Meier & Turner | [2005T ); (8) Usero et al. ( l2006l . 



0.2 0.4 0.6 0.8 




Aa(arcsec) 



Fig. 12. a) (Left panel) The Chandra X-ray image of NGC 10 68 (co ntours: 15, 30, 50, 100, 200, 400, 600, 1000, 2000 and 3000 
counts) obtained in the 0.25-7.5 keV band by Young et al. (1200 11 1 is overlaid on the PdBI SiO map (color scale in units of 
Jy km s"' beam"'), b) [Middle panel) The X-ray image obtained in the 6-8 keV band by Ogle et al. ( 2003) (contours: 0.2, 0.5, 
1, 2, 4, 8, 15, 25, 40, 80 and 100 counts) is overlaid on the SiO(2-l)/CO(l-0) brightness temperature ratio (color scale) at the 
SiO spatial resolution, c) (Right panel) Same as b) but with the X-ray image obtained in the 6-8 keV band overlaid on the CN(2- 
1)/CO(1-0) ratio (color scale) at the CO spatial resolution. Ellipses show beams of SiO and CO as in Fig.|9] 



seem to be the most successful at increasing the abundances of 
HNCO in gas phase. There is also supporting observational evi- 
dence that HNCO is related to shocks. In particular, Zinchenko et 
al. (12000) finds a good correlation between HNCO and SiO lines 
in dense cores of our Galaxy. Martin et al. ( 2008 ) have presented 
evidence of enhanced HNCO in molecular clouds suspected to 
suffer shocks in the Galactic Center. On larger scales, Meier 
& Turner (I2005I I interpret the large abundance of HNCO mea- 
sured in the nuclear spiral structure of IC 342 as due to molec- 
ular shocks. In a later study, Martin et al. (2009) observed sev- 
eral lines of HNCO in a sample of nearby starbursts and AGNs, 
concluding that molecular shocks can enhance the abundance 
of HNCO in galaxies. More recently, Rodriguez-Fernandez et 
al. (2010) have reported a significant enhancement of HNCO 
in the bipolar outflow LI 157, related to the shock chemistry at 
work in this prototypical young stellar object. 



7.1.2. Tracers of shocks in NGC 1068 

Figure [TT] shows the spectra obtained in the 3-2 and 2-1 line 
of CH3OH and in the 4-3 line of HNCO towards the CND of 
NGC 1068 with the IRAM 30m telescope. Emission is detected 
in the three lines at significant levels (S/N ratio ~7-12). Within 
the errors, the three line profiles appear centered around Wsys- 
However, the 3-2 line of methanol is a factor of 2 narrower than 
the other two lines. Such a difference cannot be attributed to 
the ~1.5x smaller beam of the CH30H(3-2) line (-17"), as the 
maximum spread of velocities is already reached on the scales 
of the CND. By contrast this result reflects the lower excitation 
of methanol lines at high-velocities (likely arising in the inner 
CND, as shown by the SiO and CN interferometer maps). This 
supports the view that there is no significant increase of molecu- 
lar densities at small radii on the scales of the CND (Sects. |6?2. II 
and 16.2.2b . 

Table 2 shows the average CH30H(2-l)/SiO(2-l), 
CH30H(3-2)/SiO(3-2), HNCO(4-3)/SiO(2-l), SiO(2- 
1)/CN(2-1) and CN(2-l)/HCN(l-0) line ratios in the CND 
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of NGC 1068 in T,„i scale. We have corrected line ratios 
from beam dilution when necessary (single-dish observations) 
using the usual prescription for a point-like source, which is 
applicable in this context for the CND. Line ratios obtained in 
NGC 1068 are compared with results obtained in the molecular 
lobes (known as Bl and B2) of the bipolar outflow LI 157 (data 
taken from Perez-Gutierrez 1999 and Nisini et al. I2007I I. We 
also list, when available, the line ratios measured in the face-on 
barred galaxy IC 342 (data taken from Hiittemeister et al. 119971 
Meier & Turner "20051 and Usero et al. 12006'). In the case of 
LI 157 we have applied a correction for beam dilution assuming 
that the sizes of molecular emission in Bl and B2 are 18" and 
30", respectively. 

As a result of this comparison we note the similarity between 
the ratios derived from CH3OH, HNCO and SiO lines in LI 157 
and NGC 1068. At face value, this indicates that most of the 
SiO emission in the CND of NGC 1068 could be explained by 
fast shocks similar to those identified in the young bipolar out- 
flow LI 157. A comparison between the SiO/CN and CN/HCN 
ratios measured in LI 157 and NGC 1068, indicates a remark- 
able excess of CN emission in NGC 1068, however. LI 157 is 
the only outflow in which CN(2-1) emission has been detected 
in the high velocity gas associated with the bow shocks. But even 
in the case of this extremely molecular-rich young outflow, the 
SiO/CN intensity ratio is ~a factor of 10 higher than in the CND 
of NGC 1068. This CN excess questions the suitability of the 
bipolar outflow template in NGC 1068. The age dating of the 
most recent star formation episode of the CND (>2-3xl0'^ yrs 
old; Davies et al. |2009l ) also questions the hypothesis that shocks 
in NGC 1068 can be interpreted as stemming from an embedded 
population of young stellar objects. 

Alternatively, large-scale shocks produced by cloud-cloud 
collisions, which are enhanced due to the complex orbital dy- 
namics in the ILR region of the CND, could be a viable mech- 
anism to explain the emission of SiO, CH3OH and HNCO in 
NGC 1068. SiO maps show the existence of perturbed kine- 
matics of molecular gas in the CND (Sect. 15.21 ). The preva- 
lence of SiO emission at high and forbidden velocities im- 
plies enhanced X(SiO) in gas following non-circular orbits. 
Furthermore, the line ratios discussed in Sect. 16.2.11 suggest an 
enhancement of SiO abundances at the extreme velocities iden- 
tified in the CND. A similar link between density waves and 
large-scale molecular shocks has been proposed to explain the 
enhancement of SiO in the circum-nuclear disks of our Galaxy 
(Martin-Pintado e tal.ll997t Huttemeister et al. 1998; Rodriguez- 
Fernandez et al. 120061) and of the starburst galaxies IC 342 
and NGC 25 3 (Use ro et al. l2006l Garcia-Burillo et al. l2000b . 
Usero et al. (2006) analyzed the large-scale molecular shocks 
produced in the ILR region of IC 342. Shocks in this galaxy 
seem to arise during cloud-cloud collisions at low velocities, 
only after the kinetic energy associated with the density wave 
driven streaming motions has partly dissipated into turbulence. 
The comparison between IC 342 and NGC 1068 indicates that 
CHaOH/SiO and HNCO/SiO ratios are ~a factor of two to three 
lower in NGC 1068. The detection of over-luminous SiO lines 
in NGC 1068 could be the signature of comparatively higher ve- 
locity shocks in this galaxy, possibly produced by mechanisms 
unrelated to density waves. Other scenarios designed to fit the 
CND kinematics, in particular those invoking a jet-ISM interac- 
tion, which would give rise to a nuclear warp and, eventually, 
to an expanding ring, could produce /oifer molecular shocks in 
NGC 1068. 

Gallimore et al. d 19961 ) analyzed the morphology and spec- 
tral index of the NGC 1068 jet and concluded that the radio 



plasma may have been diverted at a shock interface with the 
molecular ISM close to the AGN (knot C at r ~30-40pc). 
Molecular shocks in the CND could have propagated around 
the jet-ISM working surface. The H2 line map published by 
Davies et al. (2008; see also Miiller-Sanchez et al. 2009) reveals 
a strongly perturbed velocity field at a region located just 0.5- 
1" (35-70 pc) north of the AGN. The enhanced SiO abundances 
derived at high velocities 'close' to the AGN may well be the sig- 
nature of this jet-ISM interaction. Nevertheless, the fact that SiO 
gas is detected elsewhere in the CND suggests that shocks would 
be at work on larger scales, either driven by an expanding ring or 
by the perturbed kinematics at the ILR region. However, whereas 
shock chemistry can explain the strength and the line ratios of 
the different molecular shock tracers detected in NGC 1068, it is 
not a valid explanation for the measured CN abundances. 

The discrepancies mentioned above motivate the search of 
a framework different to shock chemistry that would be able to 
simultaneously account for the high molecular abundances of 
SiO and CN measured in NGC 1068 (Sect.lTIb. 



7.2. X-ray chemistry 

We describe below the main properties of X-ray emission in 
NGC 1068, with a particular emphasis on the morphology, 
energetics and origin of hard X-ray emission in the CND 
(Sect. 17.2.1b . We further discuss the observational evidence for 
X-ray chemistry in the nucleus of NGC 1068, and how the SiO 
and CN data specifically fit into this picture (Sect. l7T2!2] i. 



7.2.1. X-ray emission in NGC 1068 

We superpose in Fig. fT2h the subarcsec resolution X-ray image 
of NGC 1068 obtained in the 0.25-7.50 keV band with Chandra 
(contours) on the SiO map (color scale). The X-ray image has 
been taken from the 3.2s frame time data published by Young 
et al. (12 0011 1. The 0.25-7.50 keV band image incorporates the 
contribution from soft to hard X-rays and includes both contin- 
uum and line emission. The morphology of the X-ray image re- 
veals spatially resolved emission distributed in a compact yet 
spatially resolved component located on the AGN, and an elon- 
gated source that extends to the northeast, similarly to the inner 
radio jet of the galaxy. The northeast elongation of X-ray emis- 
sion contrasts with the overall east-west elongation seen in the 
CND molecular material detected in SiO and CN. Lower level 
X-ray emission is detected outside the inner 12"xl2" field of 
view (FOV) displayed in Fig. fTSh . This outer emission takes 
the shape of a complex spiral-like structure that extends out to 
r ~30-40" (see Young et al. 2001 for a detailed description). The 
brightest region in X-ray emission corresponds with the com- 
pact AGN source. The X-ray peak coincides with the highest 
SiO/CO and CN/CO intensity ratios measured in the CND (pan- 
els c and d). The contributions to the X-ray luminosity in the 
0.25-7.50 keV band are known to be multiple (line and contin- 
uum emission, as well as emission from neutral and ionized ma- 
terial), and therefore the link between the molecular line ratios 
and the 'bolometric' X-ray flux, if any, is not easily interpre table. 
As argued below, by restricting the study to the hard X-ray band 
the interpretation is facilitated. 

Figure. [12] (panels b and c) shows the superposition of the 
hard X-ray image obtained in the 6-8 keV band by Chandra 
(in contours) on the SiO/CO and CN/CO ratio maps (color 
scale). The X-ray map has been obtained by applying an adap- 
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Fig. 13. a) {Left panel) The SiO(2-l)/CO(l-0) brightness temperature ratio versus the Xray flux in the 6-8 keV band. The straight 
line represents the best fit to the points taking into account uncertainties visuaHzed by the errorbars. The coefficient of determination 
of the regression is p- - 0.6. b) (Right panel) Same as a) but showing the CN(2-1)/CO(1-0) ratio. The coefficient of determination 
of the regression is r^ - 0.7. 



tive smoothing □ on the Chandra archive NGC 1068 data avail- 
able in this wave-band. This image is virtually identical to the 
hard X-ray map originally published and discussed by Ogle et 
al. 



( 12003 : see also Young et al. 2001 ). Emission in this band is 
dominated by the Fe I Ka line, with a non negligible contribution 
from Fe XXV and scattered neutral and ionized continuum re- 
flection (Ogle et al. 120031) . Analysis of the 4-10 keV continuum 
observed with ASCA in NGC 1068 led Iwasawa et al. (fT997] l 
to conclude that most of the hard X-ray flux comes from reflec- 
tion by cold neutral gas. The strong emission detected in the Fe 
I Ka line in NGC 1068, a line generated by fluorescence in neu- 
tral cold molecular clouds, corroborates this picture. Therefore, 
the 6-8 keV band image allows us to directly probe to what ex- 
tent molecular gas in the CND of NGC 1068 is pervaded and 
processed by X-rays, without a significant contamination from 
other sources to the X-ray emission. 

The hard X-ray emission in the inner 12"xl2" FOV shown 
in Fig.[T2]/7,c is spatially resolved. Young et al. (201)1) confirmed 
the reality of this extended component, and concluded that the 
latter cannot be attributed to scattering by the telescope PSF. The 
compact source, located on the AGN, is suiTounded by lower- 
level emission that pervades the molecular material of the CND 
out to r ~2-3". This low-level emission shows a northeast pro- 
trusion beyond the CND extent, which aligns with the radio jet. 
The hard X-ray peak coincides within the errors with the highest 
SiO/CO and CN/CO intensity ratios measured in the CND, inter- 
preted in Sect. I6.2l in terms of a chemical enhancement (higher 
abundances) of SiO and CN. The picture of a correlation be- 
tween X(SiO), X(CN) and the degree of X-ray iiTadiation of 
molecular gas is further illustrated in Fig.[T3] There we have plot 
the SiO/CO and CN/CO ratios derived in Sect. Inversus the X- 
ray flux measured, using a 0.5" grid over the CND and a com- 
mon spatial resolution dictated in each case by the largest beam. 
Within an admittedly large scatter, Fig.[T3]provides tentative ev- 
idence of a coiTelation. The correlation is tighter in the case of 
CN: the corresponding coefficients of determination of the linear 
regressions are r^-Q.l for CN and r^=0.6 for SiO. Based on the 
larger range of molecular line ratios derived when high and low 
velocities are separately analyzed (Sect. 16.2b . we speculate that 
the correlations of Fig.[T3]would have been found to be tighter if 



using version 4. 1 of CIAO package 



the spectral resolution of X-ray observations had allowed us to 
kinematically discriminate between the two components. 

7.2.2. X-ray chemistry in NGC 1068 

The CND of NGC 1068 has the basic ingredients to become a 
giant XDR. First, the nucleus of this Seyfert 2 is a strong X-ray 
emitter. The intrinsic luminosity in the hard X-ray band has been 
estimated to be ~10'*-'-10'''*erg s"' (e.g., Iwasawa et al. 119971 
Maloney |1997t Colbert et al. 120021 ). Furthermore, the central 
engine of NGC 1068 is surrounded by a massive (~5x10^Mq) 
circum-nuclear molecular disk of ~200 pc radius. Thus, most of 
the molecular mass of the CND lies at a radius -100 pc where 
the expected hard X-ray flux is ~10-100 erg cm"^ s"'. This 
puts the molecular CND of NGC 1068 among the category of 
strongly irradiated XDRs (Meijerink et al. '2007V As argued in 
Sect. 17.2.11 the hard X-ray image of NGC 1068 provides direct 
evidence that molecular gas in the CND is pervaded by X-rays. 
We therefore expect that the chemistry of molecular gas in the 
CND should show the footprints of X-ray processing predicted 
on theoretical grounds. As discussed below, there are specific 
predictions regarding SiO and CN. 

Gas-phase models that analyze the chemistry of X-ray ir- 
radiated molecular gas foresee high abundances for CN, simi- 
lar- to those measured in the CND of NGC 1068: X(CN) ~a few 
10 **-10"^ (Lep p & Da lgai-no 1996; Meijerink & Spaans l2005l 
Meijerink et al. 120071 ). The abundance of CN is enhanced in 
XDR mainly as a result of the higher ionization degree of the 
gas. For similar reasons, there is theoretical basis supporting 
the enhancement of CN in PDR environments (e.g.; Boger & 
Sternberg 120051 Fuente et al. 1993 2008; Janssen et al. IT9951 ). 
However, based on the age estimated for the most recent star 
formation episode taking place in the CND(>2-3xlO^ yrs old; 
Davies et al. 120091 ), and considering that the UV radiation from 
young massive stars should dominate the feedback in a star- 
burst, we can safely discard PDR chemistry as the origin for the 
high abundances of CN in the CND of NGC 1068. As argued in 
Sect. 17.1.21 shocks can also be rejected as an explanation for CN 
abundances. 

The recent gas-phase models of Meijerink et al. (120071 ) 
foresee a significant enhancement of SiO in strongly irradiated 
(-10-100 erg cmr^ s"') high density (>10^cm"^) XDRs, i.e., in 
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a regime similar to the XDR environment of NGC 1068 where 
we measure X(SiO)~a few lO^'^-lO"^ . In particular, Meijerink 
et al. (2007) predict that the SiO(l-0)/CO(l-0) ratio can reach 
a value of ~0.10 assuming a standard depletion for silicon in 
gas phase. This is in close agreement with the average SiO(2- 
1)/CO(1-0) ratio observed in NGC 1068 (-0.08), assuming that 
the typical SiO(2-l)/SiO(l-0) ratios should be ~1. The latter 
is an educated guess considering the physical conditions of SiO 
gas in the CND. 

In addition to the predictions of gas-phase chemistry, 
Voit ( 119911 ) has proposed that the abundance of silicon monox- 
ide could also be boosted due to dust grain processing by X-rays. 
These can evaporate small (~10 A) silicate grains and increase 
the Si fraction in gas phase, leading to a considerable enhance- 
ment of SiO in X-ray irradiated molecular gas (Martin-Pintado 
et al. 120001 U04; Gai-cia-Burillo et al. 2008; Amo-Baladron et 
al. 120091 ). This scenario has been invoked by Martin-Pintado et 
al. C2000') and Amo-Baladron et al. (120091) to account for the cor- 
relation between the abundance of SiO and the equivalent width 
of the Fe Kc fluorescence line in the Sgr A and Sgr B molec- 
ular cloud complexes of the Galactic Center In the case of the 
Galactic Center, pure gas-phase XDR models have difficulties in 
reproducing the observed SiO abundances, unless an X-ray out- 
burst is assumed to have taken place in Sgr A* 300 yr ago. By 
contrast, current gas-phase XDR models successfully reproduce 
the SiO abundances derived in NGC 1068 without resorting to 
dust grain chemistry. Nonetheless, the inclusion of dust grain 
chemistry, likely linked to the mechanical heating of the molec- 
ular ISM, would help solve the controversy regarding the abun- 
dances of species like HCN and HNC, under-predicted by XDR 
schemes (Garcia-Burillo et al. i2008i: Loenen et al. ,2008; Perez- 
Beaupuits et al. 120091 ). Actually, while CN is overly produced in 
gas-phase XDR models, HCN is under-produced with respect to 
the level required by observations in AGNs. In addition to me- 
chanical heating, the evaporation of dust grain mantles in dense 
hot environments like those likely prevailing in the nuclear disks 
of AGN, has been invoked as a mechanism responsible of en- 
hancing HCN abundances (Lintott & Viti 2006). 

Other lines of evidence support the existence of an XDR in 
the CND. The low HCO+/HOC+ ratios measured by U04 can 
be explained only if molecular clouds have the high ionization 
degrees typical of XDR: X(e )~10"^-10 "*. Furthermore, the ex- 
citation of the 2.12yum H2 rovibrational emission Unes detected 
in the CN D (Ro taciuc et al. il991i: Bliet z et a l. IT9941 Galliano 
& Alloin I20U2I Muller-Sanchez et al. 12009] ) has been inter- 
preted to be dominated by X-ray emission (e.g.; see discussion 
in Malonev 1997 and Galhano & AUoin;2002j) 



8. Summary and Conclusions 

We have used the high spatial resolution and sensitivity capa- 
bilities of the PdBI to map the emission of the (f=0, 7=2-1) 
line of SiO and the A^=2-l transition of CN in the disk of the 
Seyfert 2 galaxy NGC 1068. The spatial resolution of these ob- 
servations (~l-3") have allowed us to separate the emission of 
the SB ring from that of the CND. The PdBI SiO and CN data, 
together with available PdBI CO maps of the galaxy, have been 
complemented by single-dish data obtained with the IRAM 30m 
telescope in SiO, CH3OH and HNCO lines to probe the physical 
and chemical properties of molecular gas in the CND through an 
analysis of line ratios. 

We summarize below the main results obtained in this work: 



- Unlike the CO lines, which show strong emission in the 
SB ring, most of the SiO emission detected inside the PdBI 
primary beam comes from a circum-nuclear molecular disk 
(CND) located around the AGN. The dichotomy between the 
CND and the SB ring is reflected in the remarkably different 
< ^sio/co > ratios measured in the two regions. LVG models 
implemented to fit line ratios indicate that the average abun- 
dance of SiO in the CND, < X(SiO)>[CND]~(l-5)xlO-'^, is 
about one to two orders of magnitude higher than that mea- 
sured in the SB ring. 

- Similarly to CO, the SiO CND, of -400 pc deconvolved 
size, shows an asymmetric double peak structure oriented 
east-west. The eastern knot corresponds to the SiO emission 
peak. However, CO and SiO emissions show differences at 
the scales of the CND. The SiO/CO velocity-integrated in- 
tensity ratio changes by about a factor of three inside the 
CND, reaching a peak ratio of ~0. 10-0. 12 towards the AGN 
and the western knot. 

- The overall SiO kinematics in the CND are consistent with 
a rotating structure. However, there is evidence that the 
CND rotating pattern is distorted by non-circular and/or non- 
coplanar motions. Compared to CO, SiO emission at small 
radii (r <2"=140 pc) stands out at extreme velocities not ac- 
counted for by circular rotation. The highest SiO/CO bright- 
ness temperature ratios in the CND are associated with the 
high velocities seen towards the AGN and the western knot. 
We interpret the spatial variations of the SiO/CO ratio in the 
CND as the likely signature of chemical differentiation. In 
this scenario, the abundance of SiO is significantly enhanced 
at high velocities out to X(SiO)~(0.3-1.5)xlO"^ 

- CN emission is also detected in a CND around the AGN. 
The size and morphology of the CND seen in CN are similar 
to those found in SiO. The overall abundance of CN in the 
CND is high: < X(CN)>[CND]~(0.2-l)xlO-^. The CN/CO 
velocity-integrated intensity ratio changes by about a factor 
of three inside the CND and reaches a peak value of ~0.30 
towards the AGN. CN maps also show distorted kinematics 
that cannot be explained by circular rotation. The abundance 
of CN is significantly enhanced at high velocities, detected 
towards the AGN, out to X(CN)~(0.8-4)xlO"'^. The emis- 
sion detected in CN at extreme velocities is physically asso- 
ciated with gas lying at small radii (r <1"=70 pc). 

- Different models have been proposed to account for the kine- 
matics of molecular gas in the CND. Large-scale shocks pro- 
duced by cloud-cloud collisions in the ILR region of the 
NGC 1068 nuclear bar can explain the distorted kinematics 
and enhanced SiO abundances in the CND. The hypothesis 
of a jet-ISM interaction, driving the large-scale expansion of 
the molecular ring, can also explain shocks. Independent ev- 
idence of shocks is provided by the detection of CH3OH and 
HNCO. Line ratios involving these tracers and SiO in the 
CND are similar to those measured in prototypical shocked 
regions in our Galaxy. However, the strength and abundance 
of CN in NGC 1068 can be explained neither by shocks nor 
by PDR chemistry. Alternatively, the high global abundances 
measured for CN and SiO, in agreement with the theoretical 
predictions of XDR models, and the conspicuous correlation 
of CN/CO and SiO/CO ratios with the irradiation of hard X- 
rays, suggest that the CND of NGC 1068 has become a giant 
X-ray dominated region. 

The extreme properties of molecular gas in the circum- 
nuclear disk of NGC 1068 result from the interplay between 
several processes directly linked to nuclear activity. Results pre- 



16 



S. Garcia-Burillo et al.: Molecular gas chemistry in AGN 



sented in this paper highlight, in particular, the footprint of 
shocks and X-ray irradiation on the properties of molecular gas 
in this Seyfert. Whereas XDR chemistry offers a simple expla- 
nation for CN and SiO in NGC 1068, the relevance of shocks 
deserves further scrutiny, however. In particular, sub-arcsecond 
spatial resolution observations, to be done at higher frequencies 
and in clear-cut molecular tracers, will be paramount to locate 
the occurrence of shocks in the CND and determine the gas 
mass and energies involved in the process. Increasing the spa- 
tial resolution at millimeter and sub-millimeter wavelengths will 
be also a key to explore in detail the correlation between molec- 
ular abundances and X-ray irradiation, tentatively identified in 
this work. Furthermore, increasing the number of molecular line 
transitions observed with high-spatial resolution will allow us 
to fit the different line ratios using multi-phase radiative transfer 
schemes. These are a better description of the multi -phase nature 
of the molecular ISM in galaxies. In particular, these models will 
properly take into account the fact that the emission of different 
lines can arise from regions characterized by very different den- 
sities and kinetic temperatures. 
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